Two-dimensional dissipative solitons supported by localized gain 
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We show that the balance between localized gain and nonlinear cubic dissipation in the two- 
dimensional nonlinear Schrodinger equation allows for existence of stable two-dimensional localized 
modes which we identify as solitons. Such modes exist only when the gain is strong enough and 
the energy flow exceeds certain threshold value. The observed solitons neither undergo diffractive 
spreading nor collapse. Above the critical value of the gain the symmetry breaking occurs and 
asymmetric dissipative solitons emerge. 



PACS numbers: 

The nonlinear Schrodinger (NLS) equation is proba- 
bly the most widely used mathematical model describing 
numerous nonlinear phenomena in the optics of Kerr me- 
dia Being supplied by gain and/or dissipation it is 
employed for modeling open systems and plays a particu- 
larly important role for description of optical patterns . 
Soliton formation described by the NLS equation and 
soliton properties depend dramatically on the dimension- 
ality of the system. In particular, in the conservative 
case, while the one-dimensional (ID) cubic NLS equa- 
tion with homogeneous coefficients has the well known 
soliton solutions, no stable localized states can exist in 
two or more dimensions Q. When homogeneous linear 
gain and nonlinear dissipation are added to the system, 
stable stationary localized structures do not exist even in 
the ID model Stabilization, however, can be achieved 
in ID case by considering a spatially localized gain 

In this Letter we demonstrate theoretically and nu- 
merically that the use of the localized gain source in the 
presence of nonlinear losses allows one to obtain stable 
dissipative solitons also in a 2D dissipative NLS equation. 
We consider the model 



(1) 



where q = q{S^, r) is a complex field, ^ is the propagation 
variable, r = (77, Q is a vector in the transverse plane, 

V = , the localized gain 7(r) is a real decaying 

function of r = |r|. It will be convenient to represent 
7 = pif{r), where pi is the gain coefficient and /(r) is a 
real decaying function with max{/(r)} = /(O) = 1. The 
positive constant a characterizes nonlinear losses. 

We notice recently increasing interest in exploring lo- 
calized gain [B^j which can be implemented in various 
physical systems, like layered structures pumped by in- 
homogeneous currents J6| , nonlinear materials with spa- 
tially localized doping 0], laser systems [1| and fluid me- 
chanics 0. AstoEq. O, it describes nonlinear response 
of semiconductor alloys where soliton formation is ob- 
servable for wavelengths below the half bandgap and the 
two-photon absorption is the dominating mechanism of 
losses at power levels required for soliton formation [lo| . 



Such materials are widely used for production of wide- 
band semiconductor optical amplifiers with high optical 
gain in the same spectral range 
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Solutions of Eq. ([T]) predicting dissipative solitons are 
searched in the form ijj = u;(r)e*''^, where h is the propa- 
gation constant, w{v) — Wr + iwi = u(r)e'^'''"-' is a com- 
plex amplitude, whose real and imaginary parts are re- 
spectively Wr and Wi, while u(r) and 9(r) are respectively 
the real amplitude and phase satisfying the equations 



- 2bu + V^u - uv^ + 2u-^ : 



0, 
= 0, 



(2a) 
(2b) 



where v = V6. We are interested in the localized solu- 
tions with zero asymptotics: u, |f | — !■ at r — > 00. 

First, concentrating on solutions with v = 0, we show 
that by controlling the gain one can construct localized 
modes of Eqs. ([2]). More specifically, we pose the problem 
of finding the gain corresponding to the a priori given 
soliton. In the ID case (r = ry and V = d/drf) for w = 
we find that Wds = e*/^^'' V [d cosh.{r] / d)] solves Eqs. ([2]) 
provided the gain is given by 7(77) = {a/ d"^) / cos\i^ {r}/ d), 
where d > is a constant defining the localization of the 
gain. Turning, to the 2D case we observe that for v = 
Eq. ([2a|) has the well known Townes soliton solution. We 
denote it as w^ir). In order to support the existence of 
the dissipative Townes soliton, the gain must have the 
form 7t(^) = aw'^{r). For any 7(r) 7^ 7t('") dissipative 
solitons with a trivial phase do not exist. 

In general, i.e. when v is not necessarily zero, station- 
ary localized solutions exist only for 5 > 0; it follows from 
the limit r ^ 00, where one can neglect both nonlinear 
terms (i.e. and u'^) in Eq. ([2a|) . The balance between 
gain and dissipation is expressed by 



^{r)\w\^ dv — a 



From ^ one can also obtain the equation 
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+^ I \Vw\'^dv = / \w\^dT 



(3) 



(4) 
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here U = J u^dr is the energy flow, that expresses the 
balance between conservative terms. Thus b is defined 
by the interplay between gain and dissipation. This im- 
poses restrictions on the range of variation of b. In- 
deed, from ([3]) and Q, combined with the estimate 
/ ^v?dY < piU, we obtain J \\7w\^dr < {pi/a — b)U. 
Hence b G (0,pi/a). 

We further take the simplest Gaussian profile 7(r) = 
Pie~^''/'^^ (other bell-shaped gain profiles yield qualita- 
tively similar results). Typical shapes of the simplest 
radially symmetric dissipative solitons obtained with the 
relaxation method are shown in Fig. [1] For a given a, 
the soliton parameters (including b) arc dictated by the 
gain coefficient. When pi is sufficiently small the solitons 
are broad [Fig. [T](a)] . The field can expand far beyond 
the spatial region where the gain is realized due to the 
energy flow outward the amplifying domain. Growth of 
Pi is accompanied by progressively increasing localization 
[c.f. panels (a) and (b) in Fig.[T]. For a > 2 with increase 
of Pi the soliton shape takes on the form of a narrow peak 
superimposed on broader beam [Fig. [T] (b)]. 
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FIG. 1: (Color online) Profiles of radially symmetric dissipa- 
tive solitons at (a) pi = 0.642, a = 1.2 and (b) pi — 4, a = 1.9. 
Real, imaginary parts of the field, and the field modulus are 
shown with black, red, and green lines, respectively. In the 
both cases d = 1.5. 

The energy flow [/ of a dissipative soliton is a non- 
monotonic function of the gain coefficient pi [Fig. [DJa)]. 
For small and moderate nonlinear losses the initial de- 
crease of U is followed by its growth, and then again 
replaced by decrease as pi increases [Fig. [2] (a)]. In the 
cutoff point p'i° the tangential line to U{p) is vertical, but 
the value If^" = U (p™) remains finite. We notice that 
it is possible to find the second branch of solutions van- 
ishing in the same point {Pi°, W") for which energy fiow 
increases with pi in the vicinity of cutoff, i.e. in (p™, W") 
two soliton families smoothly join, rather than one fam- 
ily terminates. However, this second branch is unstable 
as the linear stability analysis predicts. The character 
of U{pi) dependence changes dramatically when a is suf- 
ficiently large [Fig. [IJb)]. In this case U grows with pi 
everywhere except for a region close to the cutoff. In- 
terestingly, this change in the behavior of U{pi) accom- 
panies the tendency for development of a narrow peak 
in the center of soliton. b increases monotonically with 
Pi for small and moderate values of a, but when a ap- 



proaches 2 the dependence b{pi) becomes nonmonotonic 
[Fig.[2Ic)]. For large a one also observes two soliton fam- 
ilies merging at p^". At a fixed width of the amplifying 
domain, grows almost linearly with increase of non- 
linear losses [Fig. [2][e)]. Meantime, Pi° diminishes with 
increase of d, so that at d -H- oo we naturally recover the 
case of medium with uniform gain where solitons can be 
found for any pi [Fig. [T]^f)]. Notice however, that for 2D 
solitons considered here value also diminishes with d, 
so that while localized solitons at d — >■ oo exist for any pi 
they are all unstable. 
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FIG. 2: (Color online) Energy flow of radially symmetric dis- 
sipative soliton (a),(b), its propagation constant (c), and real 
part of perturbation growth rate (d) vs pi for different a val- 
ues and d — 1.5. Cutoff and critical values of gain coefficient 
(e) vs a at d — 1.5 and (f) vs d at a = 1.5. 

The stabilization of 2D solitons means that for the 
respective initial conditions the localized gain does not 
induce collapse, which is arrested by the nonlinear dis- 
sipation. On the other hand the nonlinear dissipation 
suppresses instability of the trivial zero solution prevent- 
ing growth of the background noise. This is in contrast 
to the case of Kerr medium with uniform gain and two- 
photon absorption where solitons are unstable even in 
ID geometries These facts have been numerically 
verified by the linear stability analysis with the ansatz 
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*g<5 5 ^ where n is an integer az- 
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imuthal index. Solution of the respective linearized equa- 
tions reveals a surprising result: most destructive pertur- 
bations correspond to n = 1, rather than to n = 0. The 
typical dependencies of real part of perturbation growth 
rate Sr on pi for n = 1 are shown in Fig. [IJd). For small 
and moderate nonlinear losses Sr vanishes for any n when 
Pi° < Pi < pTi where is some critical value. Thus 
the dissipative solitons are linearly stable in the region 
adjacent to the cutoff p™. The stability domain slowly 
expands with increase of a [Fig.^e), shaded] and broad- 
ens considerably with decrease of the width d of ampli- 
fying domain [Fig.[2{f), shaded]. At small and moderate 
a inside the instability domain Sr increases monotoni- 
cally with pi, but for sufficiently high nonlinear losses 
the character of this dependence changes [Fig. [2][d)] and 
another stability domain appears also at high pi values. 
This stability domain is limited since with further in- 
crease of Pi perturbations with n = 2 become destruc- 
tive. Results of linear stability analysis are supported by 
direct propagation of perturbed radially-symmetric dissi- 
pative solitons. Stable propagation of a relatively broad 
soliton taken from stability domain adjacent to p^" is 
shown in Fig. Ela), while stable propagation of a high- 
amplitude soliton with the shape resembling a narrow 
peak on broader pedestal from stability domain at high 
Pi and a is shown in Fig.[3]^b). 
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FIG. 3: (Color online) Stable propagation of radially- 
symmetric solitons at (a) pi = 0.8, a = 1.2, and (b) pi = 
4.2, a — 1.9, and (c) transformation of unstable radially- 
symmetric soliton into asymmetric stable state at pi — 2.0, 
a = 1.2. Field modulus distributions are shown at different 
distances. In all cases d — 1.5. 

The development of instability for perturbed unsta- 
ble radially symmetric solitons results in their transfor- 
mation into strongly asymmetric non-rotating dissipative 
soliton whose center is slightly displaced from the origin 
towards the periphery of amplifying domain [Fig. EJc)] . 
Such states can be excited with Gaussian beams too, i.e. 
they are attractors with large basin. The energy flow and 
propagation constant of the asymmetric soliton remain 
the same as long as p^, d, and a are fixed and only the 
beam orientation changes for different inputs and noise 
realizations. Thus the medium with localized gain and 



nonlinear losses can support strongly asymmetric dissi- 
pative solitons despite the fact that all parameters in the 
system are either uniform (diffraction, nonlinearity, and 
losses) or radially symmetric (gain). The asymmetric 
solitons are characterized by a specific phase distribution 
[Figs, mja) and (b)] that is asymmetric along longer axis 
of the beam and symmetric along its shorter axis. 




FIG. 4: (Color online) Field modulus (top) and phase (bot- 
tom) distributions in asymmetric dissipative solitons for pi = 
2.9 (a) and pi = 1.7 (b) at a = 1.9, d = 1.5. (c) Stable 
propagation of a perturbed asymmetric soliton at pi — 2.0, 
a = 1.2, d — 1.5. Initial and final field modulus distributions 
are shown. 

Using a method of direct propagation we obtained the 
entire family of asymmetric solitons. The dependencies of 
energy flow and integral widths = 2 (J r^^jgj^iir/ UY^'^ 
(d^ is defined similarly) on pi are shown in Fig. [5] The 
family of asymmetric solitons bifurcates from the fam- 
ily of radially symmetric solitons with increase of pi. 
The bifurcation occurs exactly in the point p"' where 
the radially symmetric solitons become unstable. With 
increase of pi the asymmetry of soliton shape increases 
[c.f. Figs, mja) and (b)]. The obtained asymmetric soli- 
tons propagate stably in the presence of small-scale input 
noise that also do not change the orientation of soliton 
in the transverse plane [Fig. UJ^c)]. 
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FIG. 5: (Color online) (a) Energy flows of asymmetric (red 
curve) and for radially symmetric (black curve) solitons vs pi 
at a = 1.9, d = 1.5. (b) drj( vs pi for the family of asymmetric 
solitons whose representative shapes are shown in Figs. |4{a) 
and (b). 
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Summarizing, wc predicted that localized gain can sup- 
port stable 2D solitons in a medium with Kerr-type non- 
linearity and strong two-photon absorption. This is in 
drastic contrast to behavior of solitons either in 2D con- 
servative systems described by the nonlinear Schrodinger 



equation, or in dissipative systems with uniform gain. 
Moreover we obtained a stable family of asymmetric non- 
rotating solitons which have no counterparts in conser- 
vative or dissipative uniform systems with cubic nonlin- 
earity. 
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